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Abstract 

In this paper we study the smah time asymptotics for the heat kernel on a sub-Riemannian 
manifold, using a perturbative approach. We explicitly compute, in the case of a 3D contact 
structure, the first two coefficients of the small time asymptotics expansion of the heat kernel 
on the diagonal, expressing them in terms of the two basic functional invariants x ^-i^d k defined 
on a 3D contact structure. 

1 Introduction 

In this paper we consider the problem of finding the small time asymptotic expansion on the 
diagonal of the heat kernel on a sub-Riemannian manifold M.| 

By a sub-Riemannian manifold we mean a triple S = (M, 2?,^), where M is a connected 
orientable smooth manifold of dimension n, 2? is a smooth vector distribution of constant rank 
k < n, satisfying the Hormander condition and g is an Euclidean structure on D. 

It is well known that a sub-Riemannian manifold has a natural structure of metric space, where 
the distance is the so called Carnot-Caratheodory distance 

d{<lo,Qi) = inf{ J \Jg'y(t){i{t),iit)) dt \ J : [0,T] M is a Lipschitz curve, 

7(0) = qo,7{T) = qi, ^{t) G a.e. in [0,r]}. 

As a consequence of the Hormander condition, this distance is always finite and continuous, and 
induces on M the original topology (Chow Theorem, see e.g. Q). 

Define inductively the vector distributions := V, := + ["D*, D], for every i = 1, 2, . . . 
and denote by T>q the fiber of at the point q. Under the hypothesis that the sub-Riemannian 
manifold is regular, i.e. for every i = l,...,m, the dimension of does not depend on the 
point, the Hormander condition guarantees that there exists a mimimal m E N, called step of the 
structure, such that = TgM, for ah q G M. 

Moreover, it is well known that the Hausdorff dimension of M, as a metric space, is given by 
the formula (see [p7| ]) 



Q = ^^ifcj, ki := dimP* — dimP 



i=l 
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In particular the HausdorfF dimension of a sub-Riemannian manifold is always bigger than its 
topological dimension. 

A sub-Riemannian structure on a three-dimensional manifold M, is said to be contact if its 
distribution P is the kernel of a differential one form P = ker w, such that A w is a nonvanishing 
3-form on M. 

The geometry of 3D contact sub-Riemannian structures have been deeply studied and they 
possess two basics differential invariants % and k (see for the precise definition and their role in 
the asymptotic expansion of the sub-Riemannian exponential map). These invariants x ^-nd k are 
smooth real functions on M that reflect geometric properties of the structure, and in the case of 
left-invariant metrics on Lie groups they have been used in ||l| to provide a complete classification 
of such structures up to local isometrics. 

Our goal is to compute the first term in the asymptotic expansion of the heat kernel in the 3D 
contact case and to understand how the invariants x k enters in the coefficients. 

The heat equation on a sub-Riemannian manifold is a natural model for the description of a 
non-isotropic diffusion process on a manifold. It is defined by the second order PDE 

d 

—Tp{t,x) = Af'il;{t,x), Vt>0, xGM, (1) 

where A j is the sub-Riemannian Laplacian (also called sub-Laplacian) which is a hypoelliptic, but 
not elliptic, second order differential operator. Locally this operator can be written in the form 
"sum of squares" plus a first order part 



A/ = ^/f + ai/j, fc = dimP<n. 



i=l 



where /i , . . . , is an orthonormal basis for the sub-Riemannian structure and ai , . . . , are suit- 
able smooth coefficients (see Section ^ for the precise meaning of these coefficients). 

From the analytical viewpoint, these operators, and their parabolic counterpart dt — Aj, have 
been widely studied, starting from the celebrated work of Hormander |24]. A probabilistic approach 



to hypoelliptic diffusion equation can be found in [10, 18, 25], where the existence of a smooth heat 
kernel for such equations is given. 

On the other hand, a "geometric" definition of the Laplacian is needed if one want to find 
some relation between the analytical properties of the heat kernel (e.g. the small time asymptotics) 
and the geometric properties of the manifold, like in the Riemannian case (see pO| , f or the 
relation between the heat kernel and the Riemannian curvature of the manifold, and [|29| for a 
characterization of the cut locus via the heat kernel) . 

As it was pointed out in ||^, |2^, to have an intrinsic definition of the sub-Laplacian Aj (i.e. 
that depends only on the geometric structure) it is necessary to build an instrinsic volume for the 
structure. 

In the sub-Riemannian case there are two intrinsic volumes that are defined, namely the Popp's 
volume (that is the analogue of the Riemannian volume form in Riemannian geometry) and the 
Hausdorff volume. In it is proved that, starting from dimension 5, in general they are not 
proportional. On the other hand in the 3D contact case they coincide. 

The existence of an asymptotic expansion for the heat kernel was proved, besides the classical 
Riemannian case, when the manifold is endowed with a time dependent Riemannian metric in 

fc(fc-l)N 



2|, in the sub-Riemannian free case (when n = k -\ — ^ — -) in ||T^. In |]26| , |l^ , |33|| the general 
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sub-Riemannian case is considered, using a probabilistic approach, obtaining different expansions 
on the cut locus and out of that. 



The same method was also applied in to obtain the asymptotic expansion on the diagonal. In 
particular it was proved that, for the sub-Riemannian heat kernel p{t, x, y) the following expansion 
holds 

p{t,x,x) ^j:^{aQ+ait + a2t^ + ... + akt^ + 0{t^+^)), for t ^ 0, (2) 

where Q denotes the Hausdorff dimension of M. 

Besides these existence results, the geometric meaning of the coefficients in the expansion on 
the diagonal (and out of that) is far from being understood, even in the simplest 3D case, where the 
heat kernel has been computed explicitly in some cases of left-invariant structures on Lie groups 
only (see [§, |ll], In analogy to the Riemannian case, one would expect that the curvature 

tensor of the manifold and its derivatives appear in these expansions. 

In this paper we compute the first terms (precisely oq and ai, referring to (||)), for every 3D 
contact structure. Our main tool is the nilpotent approximation of the sub-Riemannian structure. 
Under the regularity hypothesis, the nilpotent approximation of a sub-Riemannian structure S = 
(M, at a point q S M, is endowed itself with a left-invariant sub-Riemannian structure Sg on 
a so-called Carnot group (i.e. homogeneous nilpotent and simply connected Lie group) whose Lie 
algebra is generated by the nilpotent approximation of a basis of the Lie algebra of S, as explained 
in Section ^. 

When S is a 3D contact sub-Riemannian manifold, the nilpotent approximation Sg of S at every 
point g G M is isometric to the Heisenberg group, i.e. the sub-Riemannian structure on (where 
coordinates are denoted by g = (x,y, w)) defined by the global orthonormal frame 

h = dx + f2 = dy- ^du,. (3) 

Using this approach, we reduce the problem of computing the small time asymptotic of the 
heat kernel of the original sub-Riemannian structure to the problem of studying the heat kernel 
of a family of approximating structures at the fixed time t = 1. In such a way we present the 
original structure (and as a consequence its heat kernel) as a perturbation of the nilpotent one. 
With the perturbative method we then compute the coefficient oq, which reflects the properties of 
the Heisenberg group case (being the nilpotent approximation of every 3D contact structure), and 
the coefficient oi, where the local invariant n appears. 

The main result is stated as follows: 

Theorem 1. Let M be a 3D contact sub-Riemannian structure, with local invariants x Oind k. Let 
p{t, X, y) denotes the heat kernel of the sub-Riemannian heat equation. Then the following small 
time asymptotic expansion hold 



p{t,x,x) ^ -^{l + K{x)t + Oit^)), for t^O. 



1 

6) 

Notice that the Hausdorff dimension of a 3D contact structure is Q = 4 

2 Basic definitions 

We start recalling the definition of sub-Riemannian manifold. 



3 



Definition 2. A sub-Riemannian manifold is a triple S = (MjDjg), where 
(i) M is a connected orientable smooth manifold of dimension n > 3; 

(ii) P is a smooth distribution of constant rank k < n satisfying the Hormander condition, i.e. a 
smooth map that associates to g G M a /c-dimensional subspace T>q of TgM and we have 

span{ [Xi ,[...[X,.i,Xj]]] (q) \ X, G j G N} = T.M, V g G M, (4) 

where P denotes the set of horizontal smooth vector fields on M, i.e. 

V = {X e Vec(M) I X(g) eVg ^ q e M} . 

(Hi) gq is a Riemannian metric on which is smooth as a function of q. We denote the norm of 



a vector v £ "Dg with |?;|, i.e. |?;| 



A Lipschitz continuous curve 7 : [0, T] — )■ M is said to be horizontal (or admissible) if 

7(t) G for a.e. t G [0,T]. 

Given an horizontal curve 7 : [0,r] — )• M, the length of ^ is 

%)= ri7(i)l (5) 



JO 

The distance induced by the sub-Riemannian structure on M is the function 

diqo,qi) = inf{^(7) I 7(0) = qo,l{T) = q'1,7 horizontal}. (6) 

The hypothesis of connectedness of M and the Hormander condition guarantee the finiteness and 
the continuity of d{-,-) with respect to the topology of M (Chow-Rashevsky theorem, see for 
instance p). The function d{-,-) is called the Carnot-Caratheodory distance and gives to M the 
structure of metric space (see jT^ ). 

Locally, the pair (T>, g) can be given by assigning a set of k smooth vector fields spanning V 
and that are orthonormal for g, i.e. 

Pq = span{/i(g), . . . ,/fc(g)}, gg{fi{q),fj{q)) = 5ij. (7) 

In this case, the set {/i, • • • , /fc} is called a local orthonormal frame for the sub-Riemannian struc- 
ture. 

Definition 3. Let P be a distribution. Its flag is the sequence of distributions C T)^ C . . . 
defined through the recursive formula 

V^:=V, V'+'^ := V' + [V\V]. 

A sub-Riemannian manifold is said to be regular if for each i = 1,2,... the dimension of P* does 
not depend on the point q G M. 

A sub-Riemannian manifold is said to be nilpotent if there exists an orthonormal frame for the 
structure {/i, . . . , fk} and j G N such that [/jj, [fi^,... , [fij_j^, fij]]] = 0, for every commutator of 
length j. 
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2.1 3D contact case 



In this section we focus on three dimensional sub-Riemannian manifolds, endowed with a contact 
distribution. 

Definition 4. Let M be a smooth manifold, with dimM = 3. A sub-Riemannian structure on 
M is said to be contact if P is a contact distribution, i.e. T) = kerw, where uj G A^M satisfies 
a; A do; / 0. 

Remark 5. Notice that a contact structure is forced to be bracket generating. The fact that u is 
a contact form can be rewritten as dijj\x> 7^ 0. Moreover, it is possible to normalize the contact fom 
uj in such a way that duj\x> coincides with the volume form on V naturally defined by g (see |5|). 

Since we deal with local properties, it is not restrictive to assume that the sub-Riemannian 
structure is defined by a local orthonormal frame /i, /2 that satisfies: 

(M, w) is a 3D contact structure, 

V = span{/i, /2} = ker w, (8) 
sUi, fj) = ^ij, duj{fi, /2) = 1. 

Recall that the orthonormal frame /i,/2 is not unique. Indeed every rotated frame (where the 
angle of rotation depends smoothly on the point) defines the same structure. 

The Reeb vector field associated to the contact structure is the unique vector field /o such that 

w(/o) = 1, 

duifo, •) = 0. (9) 

Clearly, /o depends on the sub-Riemannian structure (and its orientation) only and not on the 
frame selected. 

The Lie algebra of vector fields generated by /o, /i, /2 satisfies the following commutation rela- 
tions 

[/i,/o] = cji/i + C§i/2, 

[/2,/o] = 4/i + 4/2, (10) 
[/2, fi] = c\ji + 4/2 + /o, 

where c^j are functions on the manifold, also called structure constants of the Lie algebra. 

Next we recall the definition of the local invariants of a contact three-dimensional structure. 
Here we simply give their expression when computed by means of coefficients of (|lO|). For a 
geometric interpretation of these invariants and an intrinsic definition one can see ||5|, [l| . 

Definition 6. With reference to notation (|lO|), we define 

- the first invariant 

.(.rt^OT. c.(,^,^;Jjj^,/, <=§. + f)/^), (U, 
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- the second invariant 

<1) = h{c\2) - h{cl2) - {C\2? - {Cl,f + (12) 

Remark 7. Notice that, from their very definition, x ^-iid k are smooth functions on M that do 
not depend on the choice of the orthonormal frame. Moreover x ^ and it vanishes everywhere if 
and only if the flow of the Reeb vector field /o is a flow of sub-Riemannian isometrics for M. 

3 The sub-Laplacian in a sub-Riemannian manifold 

In this section we compute the intrinsic hypoelliptic Laplacian on a regular sub-Riemannian man- 
ifold {M,T>,g), also called sub-Laplacian, writing its expression in a local orthonormal frame. In 
particular we find its explicit expression in the 3D contact sub-Riemannian case in terms of ([To|). 

The sub-Laplacian is the natural generalization of the Laplace-Beltrami operator A defined on 
a Riemannian manifold, that is A(j) = div(V</)), where V is the unique operator from C°°{M) to 
Vec(M) satisfying 

g(V(/>, X) = d(p{X), y X e Vec(M). 

Here g denotes the Riemannian metric, and the divergence of a vector field X is the unique function 
divX satisfying 

Lx/i = (divX)/i, (13) 

where ^ is the Riemannian volume form and Lx denotes the Lie derivative. 

In the sub-Riemannian case these definitions are replaced by the notions of horizontal gradient 
and of divergence with respect to the Popp measure, which is well defined in the regular case (see 

i)- 

Definition 8. Let M be a sub-Riemannian manifold and (j) £ C°^{M). The horizontal gradient of 
(f) is the unique horizontal vector field Vcp G D that satisfies 

g(V0, X) = d(t){X), y X €V. (14) 

Given a local orthonormal frame {/i, ■ ■ ■ , A} for the sub-Riemannian structure, it is easy to see 
that the horizontal gradient V(/> G P of a function is computed as follows 

k 

V</) = 5^/.(0)/., <^GC-(M), (15) 

i=l 

where the vector field acts on functions as a differential operator. 

Notation. In what follows we will denote by Aj the sub-Laplacian associated to the sub- 
Riemannian structure defined by the local orthonormal frame /i, . . . , /fc. Actually this definition 
does not depend on the choice of the orthonormal frame (see also Proposition |l^) . 

In the sub-Riemannian regular case, even if there is no scalar product defined in TgM, we can 
still define an intrinsic volume, called Popp volume, by means of the Lie bracket of the horizontal 
vector fields (see Q). Here we recall a convenient definition of Popp measure only for the 3D 
contact case. 
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Definition 9. Let M be an orientable 3D contact sub-Riemannian structure and {/i,/2} a local 
orthonormal frame. Let /o be the Reeb vector field and i'q, 1^1,1/2 the dual basis of 1-form, i.e. 
{ui, fj) = 6ij. The Popp volume on M is the forniQ fi := uq A ui A V2. 

Using the formula 

div(aX) = Xa + adivX, Va G C°°{M), X G Vec(M), (16) 
it is easy to find the expression of the sub-Laplacian with respect to any volume 

k 

div(V0) = j;div(/,(0)/,) 

1=1 

k 

1=1 

Thus 

k 

Af = Y^f! + {dwm,. (17) 
1=1 

Remark 10. Here we collect few properties of the sub-Laplacian that immediately follows from 
the definition: 

(i) The sub-Laplacian is always presented as sum of squares of the horizontal vector fields plus a 
first order horizontal part, whose coefficients heavily depends on the choice of the volume (see 
([l7| ) ) . Moreover A is the sum of squares if and only if all the vector fields of the orthonormal 
frame are divergence free. 

(ii) From (|l^) and (jl^) it easily follows that the sub-Laplacian is a homogeneous differential 
operator of degree two with respect to dilations of the metric structure. More precisely, if 
we consider the dilated structure (denoted by A/) where all vector fields of the orthonormal 
frame are multiplied by a positive constant A > 0, we have 

Axf = X^Af. (18) 



(iii) Define for every pair of functions 4',(p G (M) the bilinear form 



Jm 

Given a vector field X E Vec(M), its formal adjoint X* is the differential operator that 
satisfies the identity 

{Xc^,^)2 = {<p,x*^)2, y<p,^ec^{M). 



^In the dual basis of the frame {/i, /2, [fi, /2]} was considered to build the Popp volume. From ( p^ it follows 
that these two constructions agree each other. 
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It is easily computed that X* = —X — divX. In particular it follows that the sub-Laplacian 
is rewritten as 

k 
i=l 

and that satisfies the identities 

(A;0,99)2 = (</.,A;V9)2, ((/., A/</.)2 < 0, V ^ G Cq"" ( Af ) . (19) 

Prom this we can explicitly compute the sub-Laplacian in the 3D contact case (see also Q, 
formula (8)). 

Proposition 11. Let M he a 3D sub-Riemannian manifold and /i,/2 be an orthonormal frame 
and fo be the Reeb vector field. Then the sub-Laplacian is expressed as follows 

Af = f! + fi + cl^fi-cy2, (20) 
where 0^2 , 2 are the structure constant appearing in (|lO[) . 



Proof. Using ([LTD, it is enough to compute the functions Oi := div^/j, where fi = ui A 1/2 A is the 
Popp measure. 

To this purpose let us compute the quantity LxfJ-, for every vector field X. Recall that the 
action of the Lie derivative on a differential 1-form is defined as 



r ^ 



y, 1/ G AVM). 

i=0 



where e*"^ denotes the flow on M generated by the vector field X. Using the fact that Lx is a 
derivation we get 

Lx{i^o A z^i A 1^2) = LxVq a 1^1 a z^2 + 1^0 a Lx^i A 1^2 + t'o A i^i A LxV2- (21) 
Moreover, for every i = 0, 1, 2, we can write 

2 

i=o 

The coefficients fljj can be computed evaluating Lx^i on the dual basis 

lij = {Lxi'i, fj) 



(- 

^dt 



d tx. 



i=0 



Plugging these coefficients into (^l| ) we get 

Lx/i = {{lyi, [fl,X]) + {172, [f2,X]) + {uo, [fo, X]))fl, 



from which it foUows 

divX = {ui, [fi,X]) + {U2, [f2,X]) + {uo, [fo,X]). (22) 

Using ( p^ ) it is easy to see that {vq, [/q, X]) =0 for every horizontal vector field. Thus, applying 
(p^ with X = fi, with i = 1,2 one gets 



div/i = (z>2, [/2,/l]) = Ci2, div/2 = (z^i, [/l,/2]) 



Then (^) easily follows from (|T^). □ 

From the above construction it is clear that the sub-Laplacian depends only on the sub- 
Riemannian structure and not on the frame selected, i.e. it is invariant for rotations of the or- 
thonormal frame. Here we give also a direct proof of this fact. 

Proposition 12. The sub-Laplacian is invariant with respect to rotation of the orthonormal frame. 
Proof. Let us consider an orthonormal frame /i , /2 and the rotated one 

71 = cose/i + sin0/2, (23) 

72 = -sin^/i + cos^/2, (24) 

where 9 = 9{q) is a smooth function on M. If we denote by c^^ the structure constants computed 
in the rotated frame, from the formula 

[I1J2] = [/i,/2] - fi{e)fi - f2{e)f2, (25) 

it is easy to prove that the new structure constant are computed according to the formulas 

c\2 = COSe{c\2 - fl{d)) - Sin^(c?2 - f2{e)), 

= sin0(cl2 - fi{9)) + cos^(c?2 - f2{e)). 
From these relations one can easily compute 

/i + 71 = /i + /I + fmf2 - /2(^)/i, (26) 

and 

- C\2J2 + 2?27 = -(cL - + (C?2 - /2(^))/l. (27) 

Combining ( p6| ) and (27) one gets, denoting Aj the laplacian defined by the rotated frame, 

Aj=7l +7|-cl272 +2'?27 

= fl + /I + fmf2 - f2{e)fl - cl^h + 2?27 

= f! + /I + fli0)f2 - f2{e)fl - {c\2 - fl{e))f2 + (C?2 - /2(0))/l 
= /l +/|-cl2/2 + 4/l = A/. 

□ 
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Notice from (^) that the sum of squares is not an intrinsic operator of the sub-Riemannian 
structure. 

Remark 13. The same argument provides a proof of the fact that, on a 2-dimensional Riemannian 
manifold M with local orthonormal frame /i , /2 that satisfies 

[/i,/2] = ai/i + 02/2, 

the Laplace-Beltrami operator is locally expressed as 

A = /f + /I + 01/2-02/1. 

4 The nilpotent approximation 

In this section we briefly recall the concept of nilpotent approximation. For details one can see 



4.1 Privileged coordinates 

Let S = {M,T>,g) be a sub-Riemannian manifold and (/i,...,/^) an orthonormal frame. Fix 
a point q £ M and consider the flag of the distribution C 2?^ C . . . C f^. Recall that 
ki = dimP* — dimD*~^ for i = 1, . . . ,m, and that ki + . . . + km = n. 

Let Og be an open neighborhood of the point q £ M. We say that a system of coordinates 
ip : Og ^ M"" is linearly adapted to the flag if, in these coordinates, we have ip{q) =0 and 

Mvi) = r''' (B ...em'"', Vi = i,...,m. 

Consider now the splitting M" = © . . . ©M'^™ and denote its elements x = {xi, . . . , Xm) where 
Xi = (xj, . . . , S M'^'. The space of all differential operators in with smooth coefficients forms 
an associative algebra with composition of operators as multiplication. The differential operators 



J a 

a vviuii gciicicxuwio j.,X 

i = 1, . . . ,m; j = 1, . . . ,ki. We define weights of generators as 



with polynomial coefficients form a subalgebra of this algebra with generators l,x-,^~, where 

* ox-' 



= 0, u{xl) = i, 

ox. 



and the weight of monomials 



/9 



p i=l j=l 

Notice that a polynomial differential operator homogeneous with respect to v (i.e. whose monomials 
are all of same weight) is homogeneous with respect to dilations 5t : M" — )• M"' defined by 

6t{xi,...,Xm) = {txi,fx2,...,t'^Xm), t>0. (28) 

In particular for a homogeneous vector field X of weight h it holds 

St^X = t-^X. (29) 
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A smooth vector field X G Vec(M"'), as a first order differential operator, can be written as 

^ d 



and considering its Taylor expansion at the origin we can write the formal expansion 



h=—m 

where X*^'^-* is the homogeneous part of degree h oi X (notice that every monomial of a first order 
differential operator has weight not smaller than — m). Define the filtration of Vec(M") 

= {X e Vec(]R") : X« = 0,Vi < /i}, £ E Z. 

Definition 14. A system of coordinates ip : Oq ^ M" defined near the point q is said privileged for 
a sub-Riemannian structure S if these coordinates are linearly adapted to the flag and such that 
G 'D^~^^ for every i = 1, . . . ,k. 



Existence of privileged coordinates is proved, e.g. in [0, 14 1. Notice however that privileged 
coordinates are not unique. 

Definition 15. Let S = (M, P,g) be a regular sub-Riemannian manifold and {/i, . . . , fk} a local 
orthonormal frame near a point q. Fixed a system of privileged coordinates, we define the nilpotent 
approximation ofS near q, denoted by Sg, the sub-Riemannian structure on M" having {/i, . . . , fk} 
as an orthonormal frame, where fi := {ip*fi)^~^^- 

In what follows we often omit the coordinate map in the notation above, by writing fi := f^ 

Remark 16. It is well known that under the regularity hypothesis, Sq is naturally endowed with 
a Lie group structure whose Lie algebra is generated by left-invariant vector fields fi, ■ ■ ■ , fk- 
Moreover the sub-Riemannian distance d in Sq is homogeneous with respect to dilations 6t, i.e. 
d{5t{x),6t{y)) = td{x,y). In particular, if Bq{r) denotes the ball of radius r in S^, this implies 

StiBqil)) = Bq{t). 

The following Lemma shows in which sense the nilpotent approximation is the first order ap- 
proximation of the sub-Riemannian structure. 

Lemma 17. Let M be a sub-Riemannian manifold and X G Vec(M). Fixed a system of privileged 
coordinates, we define X^ := edi^X. Then 

X'^ = X + eY^ , where is smooth w.r.t. e. 

Proof. Since we work in a system of priviliged coordinates, in the homogeneous expansion of X 
only terms of order > 1 appear. Hence we can write 
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Applying the dilation and using property ( p9| ) we get 

£ fc 

Multiplying by e and using that, by definition, X = X^"^), we have 

+ + e^X^^) + . . . (30) 

□ 

In other words the nilpotent approximation of a vector field at a point q is the first meaningful 
term that appears in the expansion when one consider the blow up coordinates near the point q, 
with rescaled distances. 

Remark 18. If S = (Af, is a 3D contact sub-Riemannian manifold, then dimPg = 2 and 
dimPg = 3 for all q G M. Under this assumption the nilpotent approximation Sq of S at every 
point g E M is isometric to the Heisenberg group, since this is the only nilpotent left-invariant 
structure with ^(S) = (2,3) (see e.g. for a classification of left-invariant structures on 3D Lie 
groups). 

The sub-Riemannian structure on the Heisenberg group is defined by the global orthonormal 
frame on M'^ (where coordinates are denoted hy q = {x, y, w)) 

fi =dx + |«9«„ h = dy- |(9^. (31) 

Notice that the Lie algebra Lie{/i,/2} is nilpotent since [/i, [/i,/2]] = [72, [/15/2]] = 0. More- 
over the Reeb vector field is /o = and the local invariants of the structure are identically zero 
X = = 0. 



4.2 Normal coordinates 

In the 3D contact case there exists a smooth normal form of the sub-Riemannian structure (i.e. 
of its orthonormal frame) which is the analogue of normal coordinates in Riemannian geometry. 
This normal form is crucial for the study of the heat kernel of the sub-Laplacian with a pertur- 
bative approach, since it presents the sub-Riemannian structure of a general 3D contact case as a 
perturbation of the Heisenberg (nilpotent) case. 

Theorem 19 (|^ ^). Let M he a 3D contact sub-Riemannian manifold and /i, /2 a local orthonor- 
mal frame. There exists a smooth coordinate system (x, y, w) such that 



fi = {dx + |<9^) + f3y{ydx - xdy) + ^yd^, 



X 

/2 = {dy - -dw) - /3x{ydx - xdy) + jxdw, 

where /3 = f3{x,y,w) and 7 = j{x,y,w) are smooth functions that satisfy the following boundary 
conditions 

mO,w) = 7(0,0,7/;) = 1^(0,0,//;) = |^(0,0,w;) = 0. 
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Notice that the normal coordinate system is privileged at 0. Indeed from the explicit expression 
of the frame it immediately follows that these coordinates are linearly adapted at since 

Vo = span{/i(0), ^(O)} = span{c>,, dy} = M^. 

Moreover the weights of the coordinates (x, y, w) at the origin are 

iy{x) = u{y) = 1, u{w) = 2, 

and every homogeneous term of the vector fields /i, /2 has degree > — 1. 

Finally, notice also that, when /3 = 7 = 0, we recover the Heisenberg group structure (|3l|). 



5 Perturbative method 

In this section we consider the sub-Riemannian heat equation associated to a sub-Riemannian 
structure / on a complete sub-Riemannian manifold M, i.e. the initial value problem 

-^{t,x) = Afi;{t,x), in (0,00) xM, 

(32) 

^(0, x) = ip{x), x£M, ipe C^{M). 

where "0(0, x) = Wmij^it^x) and the limit is meant in the distributional sense. 

Recall that a differential operator C is said to be hypoelliptic on a subset ^7 C M if every 
distributional solution to Cu = (p \s C°^{U), whenever cf) G C°°{U). The following well-known 
Hormander Theorem gives a sufficient condition for the hypoellipticity of a second order differential 
operator. 

Theorem 20 (Hormander, ||2^). Let C be a differential operator on a manifold M , that locally in 
a neighborhood U is written as 

k 
i=l 

where Xq,Xi, . . . , E Vec(Af). // LiegjXo, Xi, . . . , Xk} = TgM for all q £ U, then C is hypoel- 
liptic. 

From this Theorem and the bracket generating condition it follows that the sub-laplacian is 
hypoelliptic. Moreover, since M is complete and the sub-Laplacian is symmetric and negative with 
respect to the Popp's measure (see (|^)), it follows that Aj is essentially self-adjoint on C^{M) 
(see also (iii) in Remark 10). 

As a consequence the operator A/- admits a unique self-adjoint extension on L^(Af) and the 
heat semigroup {e^^f}t>o is a well-defined one parametric family of bounded operators on L'^{M). 
Moreover the heat semigroup is contractive on L'^{M) (see |3^). 

The problem (32) has a unique solution, for every initial datum Lp G L'^{M), namely Tp{t,x) := 



e^'^fip. Due to the hypoellipticity of Aj, the function (t,x) 1— )■ e^'^fip{x) is smooth on (0, 00) x M 
and 

e*^/(/;(x) = / p{t,x,y)^{y)dy, ^ G C^{M), 
where p{t, x, y) is the so-called heat kernel associated to e^'^f , that satisfies the following properties 
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(i) p{t,x,y) G C7~(M+ x M x M), 

(ii) p{t,x,y) =p{t,y,x), Vt > 0,Vx,y G M, 

(iii) x, y) > 0, Vi>0,Vx,yEM. 

A probabilistic approach to hypoelliptic diffusion equation can be found in jiol , |l8|, |2^, wliere the 
existence of a smooth heat kernel for such equations is given. 

In particular, since e^^fip satisfies the initial condition, it holds 



/ P(t,x,y)ip{y)dy = ip{:^ 
Jm 



lim 

and the heat kernel p{t, x, y) is a solution of the problem ( |32D with the initial condition ^/'(O, x) = 
Sy{x), where 6y denotes the Dirac delta function. 

For a more detailed discussion on the analytical properties of the sub-Riemannian heat equation 
and its heat kernel one can see |31, 12 1. 

To study the asymptotics of the heat kernel associated to the sub-Riemannian structure defined 
by / near a point q £ M, we consider the approximation of the sub-Riemannian structure (cfr. also 
Lemma 17). 

Definition 21. Let fi, ■ ■ ■ , fk be an orthonormal frame for a sub-Riemannian structure on M and 
fix a system of privileged coordinates around the point q G M. The e-approximated system at q is 
the sub-Riemannian structure induced by the orthonormal frame /f , . . . , /| defined by 

fi := e5i/^^fi, i = l,...,k. 

Remark 22. Notice that in the definition of approximated structure we have to perturbate the 
basis of the distribution only. Their Lie brackets are changed accordingly to the formula 

In particular, in the 3D contact case, the Reeb vector field /g of the e-approximated structure is 
related to the unperturbed one by /q = e^5i/e*/o- 

The following Lemma shows the relation between the heat kernel of / and the one defined by 
whose existence is guaranteed by the above results. 

Lemma 23. Let M be a sub-Riemannian manifold and fix a set of privileged coordinates in a 
neighborhood N of q. Denote by f the sub-Riemannian structure and by f^ its e- approximation at 
q. If we denote by p(t,x,y) and p^{t,x,y) the heat kernels respectively of the equations 

d'ip dtp 

— (t, x) = Afil;{t, x), -^{t, x) = Af.ipit, x), 

we have that 

p'{t,x,y) = e^p{e^t,5ex,6sy), yx,y£N, 
where Q denotes the Hausdorff dimension of M. 
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Proof. Recall that, if fi,---,fk is an orthonormal frame for the sub-Riemannian structure, the 
approximated system is defined by vector fields 

fi-=^h/e*fu i = l,...,k. 

Then the proof follows from the following facts 

(i). If we perform the change of coordinates x' = (Ji/gX and we denote by q^{t,x,y) the heat 
kernel written in the new coordinate system (which depends on e), we have the equality 

Indeed it is easy to see that q'^{t,x,y) is a solution of the equation written in the new variables 
x' = 6i/sX and the factor e*^ comes from the fact that | det (J^j = . More precisely if we set 

'ip%t,x)= q^t,x,y)ip''{y)dy, 
Jm 

where (p^{x) = ip{5sx), one can see that the initial condition is satisfied: 
lini / q^{t,x,y)ip^{y)dy = liin e'^p{t,6eX,5ey)ip{6sy)dy 

lim / p{t,6eX, z)(p{z)dz {z = 5ey) 



M 



= ip{6£x) = ip^x). 
{ii). Since 5i^fi = \ff, for every i = l,... ,k, using (||) we get 

This equality can be rewritten as 

and performing the change of variable t = e^r we get 

from which it follows 

p%t, X, y) = q^e'^t, x, y) = e^p{e^t, ^x, 5ey). 



□ 



This result is useful for the study of the asymptotics on the diagonal since, by construction, the 
intial point is fixed for the dilation Se . In particular we can recover the small time behaviour of the 
original heat kernel from the approximated one 

Corollary 24. Following the notations introduced above 

p^(l,0,0) =£^(^2^0,0). 
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As a corollary we recover an homogeneity property of the heat kernel for a nilpotent structure. 

Corollary 25. Assume that the sub-Riemannian structure is regular and nilpotent. Then the heat 
kernel p{t, x, y) of the heat equation satisfies the homogeneity property 

p{t, X, y) = \^p{\h, 6xx, 6xy), V A > 0. 



This result was already well-known (see e.g. In our notations, it is a consequence of the 

fact that, if / = /, then / = /*" = / for every e > 0. 

5.1 General method 

In this section we briefly recall the perturbative method for the heat equation presented in Chapter 
3 of |3C]. For further discussions one can see also 20 1 and [34]. 



Let X, Y be operators on a Hilbert space of functions. Moreover assume that X and X + Y 
(where Y is treated as a perturbation of X) have well defined heat operators e''-^ , ^ i.e. a 

semigroup of one parameter family of bounded self-adjoint operators satisfying 

(dt - X)e^^ip = 0, lim e^^ ip = 

t-)-0 

and similarly for X + Y. Given A{t),B{t) two operators on the Hilbert space, if we denote their 
convolution as 



{A*B){t) = [ A{t- s)B{s)ds, 
Jo 



then the classical Duhamel formula 



AX+Y) ^ ^tX ^ f\it-snx+Y)y^sX^^^ 

Jo 



can be rewritten as follows 

Iterating this construction one gets the expansion 

^t(x+Y) ^ ^tx ^ ^tx ^ y^tx ^ ^t(x+Y) ^ ^Ye'^'r^ (33) 

where A*"^ = A* A denotes the iterated convolution product. 

If A{t) and B{t) have heat kernels a(t, x, y) and 6(t, x, y) respectively, then [A*B){t) has kernel 
(see again |0|, ||]) 

[a * b)(t,x,y) := / / a{s,x, z)b{t — s, z,y)dzds, 

Jo J M 

Interpreting (^) at the level of kernels, denoting by pit, x, y) the heat kernel for the operator X 
and by p^ {t, x, y) the kernel of the perturbed operator X -|- y we can write the expansion 

{t,x,y) = p{t,x,y) + {p*Yp){t,x,y) + {p^ *Yp*Yp){t,x, y) (34) 
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6 Proof of Theorem J] 

In this section we compute the first terms of the small time asymptotics of the heat kernel. To this 
extent we compute the sub-Laplacian associated to the approximated sub-Riemannian structure 
and we use the pertubative method of Section 5A to compute this terms using the explicit expression 
of the heat kernel in the Heisenberg group. 

Remark 26. The sub-Laplacian on the Heisenberg group is written as the sum of squares (cfr. 
also Remarks IC and 18) 

Aj = + ^ = a - ^d^f + {dy + ^d^f. 

The heat kernel for Aj has been computed explicitly for the first time in [^3|. Here we use the 
expression given m g in the same coordinate set. Denote by g = {x,y,w) € a point in the 
Heisenberg group. The heat kernel H{t,q,q'), is presented as 

H{t,q,q') = ht{q' oq~^), (35) 

where 

and o denotes the group law in H3 

(x, y, w) o {x', y', w') = {x + x',y + y' , w + w' + ^{x'y - xy')). 
Notice that the inverse of an element with respect to o is 

{x,y,wy^ = {-x,-y,-w). 
For a discussion on the convergence of the integral ( |36| ) one can see [^]. 

6.1 Local invariants 

In this section we compute the invariants x *iid k at the origin of the sub-Riemannian manifold. 
By Theorem 19 we can assume that the orthonormal frame has the form 

fi = {dx + l^^) + /3y(.ydx - xdy) + 7^(9^, 

/2 = {dy - -dy,) - Px{ydx - xdy) + jxd^, (37) 
where /3 and 7 are smooth functions near (0, 0, 0) that satisfy 

m 0, w) = 7(0, 0, w) = 1^(0, 0, w) = ^(0, 0, w) = 0. (38) 



Moreover, since we are interested up to second order terms in the expansion (40), we can assume 
the following 
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Lemma 27. We can assume that the orthonormal frame has the form 

/i = a,-|(i + 7)5^, 

/2 = 5, + |(l + 7R, (39) 

where j is a quadratic polynomial of the form "/{x, y) = ax^ + bxy + cy^, for some a,b,c ^ M. 

Proof. Recall that if we expand a vector field X in homogeneous components (when written in a 
privileged coordinate system) 

its e-approximation has the following expansion 

~ X + eX(°) + e^X^i) + e^X^^) + . . . (40) 

It is then sufficient to consider, in the Taylor expansion of the orthonormal frame near the 
origin, only the homogeneous term up to weight one, since every other term with weight > 2 gives 
a contribution o(e^) when one compute the expansion of Aje. Hence a contribution o(e^) in the 
heat kernel due to (^). 

Moreover the boundary condition (|3^ ) implies that the following derivatives of the coefficients 
of (^) vanish at the origin 

1^(0,0,0) = g(0,0.0) = ^(0,0,0) = ^(0.0,0) =0, (41) 

together with all higher order derivatives with respect to w. 

Since ^{dx) = v{dy) = 1 and /3(0, 0,0) = 0, the terms j3y{ydx — xdy) and f5x{ydx — xdy) have 
weight > 2. Moreover y{dw) = —2 implies that xdyj and ydyj have weight —1. The only terms that 
we need in the expansion of 7 are those of weight less or equal than one. Since the terms of order 
zero vanish by (]3^ ) and (|4l|) , the only meaningful term in the expansion of 7 is 



, , (9^7 2 d"^! 2 

7(x, y, w) ~ ^-^x +- — xy + —ry 
ox'^ oxy oy^ 



where derivatives are computed at the origin (0,0,0). □ 
Now we express the invariants x &iid k in terms of the perturbation (^9|). 



Lemma 28. Assume that the orthonormal frame of the sub-Riemannian structure has the form 
Then value of the invariants at the origin are 



X 



2v^62 + (c- a)2, K = 2{a + c). (42) 



Proof. To compute the invariants we need to compute the Reeb vector field /o and the structure 
constant of the Lie algebra Lie{/o, /i, /2}- Every contact form for the structure is a multiple of 

uj = dz - -{I + -f)dy + 2 (1 + 'y)^^^ 
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whose differential is computed as follows 

dcD = — (1 + 2^)dxdy. 

Since (iw(/i,/2) = ~(1 + 27), the normalized contact form oj that satisfies duj{fi,f2) = 1 (see 
Remark and (0) ) is 



1 



CO :- 



1+27 

Notice that every contact form vanishes on the distribution. Thus 
Next we compute the differential of the normalized contact form 



d{l + 2j)Auj , , 
du; = — ^^7: h dxdy 



(1 + 27)2 
29x7 



(1 + 27)2 



dxdw + 



2dyj 
(1 + 27)2 



dydw + ( 1 



27(1+7) 
(1 + 27)2 



dxdy, 



The Reeb vector field is, by definition, the kernel of duj (normalized in such a way that ^(/o) = 1) 
From this one gets 



/o 



2^x7 ^ _ 29,7 , 
1 + 27"^^ 1 + 27"^" 



27(1 + 7) 



(1 + 27) 9, 



1 + 27 

The commutator between horizontal vector fields is computed as follows 



[/2, fi] =[dy + ^{l + 7)9^,9.. - |(1 + 7)9„ 



= -(1 + 27)9^, 

and writing [/2, fi] = fo + ^12/1 + C12/2 we find the structure constants 

2dyj 2 29a;7 



C12 



^12 



1 + 27' 1 + 27' 

Moreover, a longer computation for [/i,/o] and [/2,/o] shows that 

2 



(1 + 27) 



r' 



2 Oi 



where we set 



Cqi = (1 + 2j)dxy-/ - 2dyjdxj, 
-(l + 27)a,,7 + 2(9,7)2, 



7^ 



2-^2 = (1 + 27)93^,7 - 2(9,7)' 



?2 
''02 



-(1 + 27)9^,7 + 29,79,7. 



19 



Recalling that at the origin dx"f = dy^ = 0, while dxxl = 2a, dxy^ = b, dyy^ = 2c it follows from 



(|n|) that 



X = 2y-det 

and 



h c — a 
c — a —b 



) = 2,/b^ + {c-ay. 



^ = /2(cL) - /i(c?2) - icl2? - {cl,f + '-^^ = 2(a + c). 



6.2 Asymptotics 

In this section we compute the Laplacian Aje up to second order in e. First notice that 

ri = {d.-ld^)-e\lld^) + o{e'), 



Moreover, defining /q := 6^5i^/o (see Remark 22), from the formula 

[/i,/f] = /o^ + (cLr/f + (c?2r/i, 

we get the following expansion 



□ 



{cl,r = e^^^ = 2e%j + o{e% (43) 



(cl^r = -^f^ = -2^'^^7 + oie'). (44) 
Thus we can compute every term defining the sub-Laplacian 

if!? = (dx - ^d^f - {{dx - |5^)(|79») + i^ld^W. - |a^)) + o{e^) 

2 

= (fif - e^{y-fdy,x - Y^dl + ^dxjd^) + o(e^), 
ifl? = idy + |9»)2 + [{dy + |5.)(|75.) + (|79.)(5, + |9.)) + o{e') 

X 

= {hf + e^{x-idyjx + -i^lSli + -dy-fdu,) + o(e^), 
where /i,/2 denote the orthonormal frame of the Heisenberg group (see again Remark 18). More- 



over, from (|4^) and (^) one easily gets 



{cl,rf! = -2e'dxi{dx - |a^) + o(e^). 
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Recollecting all the terms we find that 
where 

= + it = (5. - + {dy + 

is the sub-Laplacian on the Heisenberg group, and 3^ denotes the second order differential operator 
y = ■^{^'^ + y^)9i + lixd^y - ydyjx) - ^{xdyj - ydx7)dw - 2{dxjdx + dyjdy). 



Specifying expansion (|34| ) to our case, where Y = Y{e) is a smooth perturbation which expands 
with respect to e as follows 

Y{e) = £^y + o{£^) 

we find that 

= H + e^{H*yH) + e^{p'*yH*yH). (45) 

Since p*" is the heat kernel of a contraction semigroup (for every e > 0), one can see that the 
last term in ( ^5|) is bounded and 

p'{l,0, 0) = H{1, 0, 0) + e^H * yH{l, 0, 0) + 0{e^), 

where 

yH{t,q,q')=y,H{t,q,q'), 

means that y acts as a differential operator on the first spatial variable. 
From the explicit expression (^) it immediately follows that 

H{1,0,0) ^ 



16*2' 



Thus, denoting by Ki := H*yH{l, 0, 0) from Corollary ^ we have the expansion of the original 
heat kernel 

p{t,x,x) ~ -^{l + K,t + 0{t^)). 
We are left to computation of the convolution between H and YH, namely 

H *yH{l,0,0) = [ f H{s,0,q)yH{l- s,q,0)dqds 
Jo Jr3 

/ hs{q)yhi-s{q)dqds, q={x,y,w). (46) 
Computing derivatives under the integral sign one gets 

dr, 



(47rt)2 sinhr \ 2ttanhr / t"^ 
7(x,y)cos(— ) lr{x +y )- ^^^^^^ 1+^7 {x,y)sm.[—) 
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where 

7(3;, y) = ax^ + bxy + cy^, l'{x^ y) = 2(a — c)xy + b{y^ — x^). 

Notice that we are interested only in computing how the integral ( |46| ) depends on the constants 
a,b,c, and the perturbation of the metric is on the variables x,y only. Hence, using that the 
integrand has exponential decay with respect to x, y, we can exchange the order of integration in 
and integrate first with respect to these variables. 
Using that, for a > 

x'e"'^(^'+y'Uxdy = ^, [[ x^e-"(-^+^^) = ^, 
a JJk.2 a-^ 

xye-'^^^'+y'Uxdy = 0, 

for some constants ci, C2 > 0, it is easily seen that integrating ( ^6|) we get an expression of the kind 

H*yH{l,0,0) = Coia + c) = Cofi, (47) 

where Co is a universal constant that does not depend on the sub-Riemannian structure. Hence 
the value of Cq can be computed from some explicit formula of the heat kernel in the non nilpotent 
case. Using the expression given in |jll| for the heat kernel on SU{2), where the value of local 
invariants are constant x = 0, k = 1 (see also [||) we get that 

Psui2){t, 0' 0) = ~ 1^2 (1 + * + 0{t')), 

where we renormalized the constants in order to fit into our setting. Hence Co = 1 in (|47|), and the 
Theorem is proved. 

Remark 29. The same method applies to get a quick proof of the following- well known result (see 



e.g. |17, ^]): on a 2-dimensional Riemannian manifold M, the heat kernel p{t,x,y) satisfies an 



asymptotic expansion on the diagonal 

p{t,x,x)^-^(l + ^P-t + 0{t'')), fort^O, 

where K{x) denotes the gaussian curvature at the point x G M. Indeed one can use the normal 
coordinates on M to write the orthonormal frame in the following way 

fi = dx + f3y{ydx - xdy), 
f2 = dy- f3x{ydx - xdy), 

where (3 is, a priori, a smooth function /3 = f3{x,y). Reasoning as in Lemma /3 can be chosen 
as a constant since we are interested only in first order term. In this case it is also easily seen that 
the Gaussian curvature at the origin is computed via the parameter /3 as K = 6(3. 

Aknowledgements. The author is grateful to Andrei Agrachev for suggesting the problem 
and many helpful discussions. 
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